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Abstract. We examine the minimal magnitude of perturbations necessary to 
change the number N of static equihbrium points of a convex solid K. We 
call the normalized volume of the minimally necessary truncation robustness 
and we seek shapes with maximal robustness for fixed values of A'^. While the 
upward robustness (referring to the increase of N) of smooth, homogeneous 
convex solids is known to be zero, little is known about their downward ro- 
bustness. The difficulty of the latter problem is related to the coupling (via 
integrals) between the geometry of the hull bd K and the location of the center 
of gravity G. Here we first investigate two simpler, decoupled problems by ex- 
amining truncations of bd K with G fixed, and displacements of G with bd K 
fixed, leading to the concept of external and internal robustness, respectively. 
In dimension 2, we find that for any fixed number TV = 25, the convex solids 
with both maximal external and maximal internal robustness are regular 5- 
gons. Based on this result we conjecture that regular polygons have maximal 
downward robustness also in the original, coupled problem. We also show 
that in the decoupled problems, 3-dimensional regular polyhedra have maxi- 
mal internal robustness, however, only under additional constraints. Finally, 
we prove results for the full problem in case of 3 dimensional solids. These 
results appear to explain why monostatic pebbles (with either one stable, or 
one unstable point of equilibrium) are found so rarely in Nature. 



1. Introduction 

Ever since the work of Archimedes [I , the study of equihbrium points of a convex 
sohd, with respect to its center of gravity G has been a fundamental question of 
statics. The number N of equihbria is characteristic of the shape, it has been 
apphed to classify turtle shells [T^ and beach pebbles [12] • These classifications 
are based on the number and type of equilibrium points of a convex, homogeneous 
solid; in 2 dimensions we have S stable and S — U unstable points, and it was 
shown [11 that any equilibrium class {S} is non-empty if S" > 1. In 3 dimensions 
we have, in addition to the previous two types of equilibria, H saddle points; based 
on the Poincare-Hopf Theorem [5] we have S + U — H = 2. It is known [53] that 
equilibrium classes {5*, U} are non-empty for any S,U > 0. It is a natural question 
to ask how difficult it is to change the equilibrium class of a convex solid; we call 
this property robustness, and our aim is to introduce possible approaches to this 
concept. The notion of robustness is physically motivated by erosion processes [3], 
where small amounts of material are being abraded by collisions [151 , [5^ , [T7] and 
friction [8j. 
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The full mathematical problem can be defined by asking for the minimal nor- 
malized volume of a truncation of a convex, homogeneous solid under which N 
will change; we call this scaled volume the robustness of K and denote it by p{K). 
(To avoid confusion, we sometimes will refer to this quantity as the full robustness.) 
Depending on the sign of the change, we may define upward and downward robust- 
ness. A recent result |10l shows that the upward robustness of generic, smooth 
convex bodies (i.e. convex Morse functions on the sphere, cf. [iH]) is zero. This 
phenomenon is related to the existence of additional numbers of equilibria on finely 
discretized curves and surfaces [9] and is also distantly related to Zamfirescu's result 
|25| on the existence of infinitely many equilibria on the boundary of a typical (and 
thus neither smooth nor polyhedral) convex body. Nevertheless, the downward 
robustness (i.e. the relative volume of the smallest trunction of a smooth convex 
body decreasing N) is apparently not zero. We call the supremum of the downward 
robustness in any equilibrium class the robustness of the equilibrium class and de- 
note it by PS and ps^u in 2 and 3 dimensions, respectively. Finding the robustness 
of an arbitrary equilibrium class appears to be a nontrivial problem. Beyond being 
mathematically challenging, it also looks rather interesting from the point of view 
of natural abrasion processes. Our goal in this paper is to deliver some partial 
results which may serve as a basis for the intuition about the solution of the full 
problem. 

The main difficulty of the latter lies in the nontrivial coupling via integrals be- 
tween the hull bd K and the center of gravity G. In the first part of this paper we 
solve two simpler, decoupled problems. In the first case we seek internal robust- 
ness, defined as the minimal (normalized) distance necessary to move G at fixed 
hdK leading to a change in N. The concept of internal robustness is phyiscally 
motivated by material inhomogeneities. From the mathematical point of view, this 
problem is closely related to the geometry of caustics 19!, and has attracted re- 
cent interest in the context of inhomogeneous polyhedra (cf. and [H]). Here 
we do not distinguish between upward/downward robustness, rather we compute 
their minimum. The other decoupled problem leads to the concept of external ro- 
bustness, defined as the minimal (scaled) truncation of K at fixed G, leading to 
a change in N. Similarly to the full problem, (in fact, as a consequence of the 
same theorem [TD]) here again the upward robustness is zero. We show that in the 
plane, both decoupled problems lead to the same result: for any fixed N = 2S, the 
convex shapes with maximal internal or external (downward) robustness are reg- 
ular S-gons, thus we determine the internal and external robustness for all planar 
equilibrium classes. This result suggests that regular polygons may have maximal 
robustness in the original, full problem as well. In the case of internal robustness, 
we show also that platonic solids have maximal robustness, however, only under an 
additional constraint. After exploring the decoupled problems we investigate the 
full (downward) robustness of some selected equilibrium classes in 3 dimensions, in 
particular, we show that if S,U < 3 then ps,c/ = 1- We also explore the partial 
robustness of convex 3D solids, measuring the difficulty to either reduce S or U. 
Our results offer one plausible explanation to the geological puzzle, why monostatic 
pebbles (with cither S = 1 or U = 1) are found so extremely rarely in Nature. 

After introducing basic notions and notations in Section [2j we investigate exter- 
nal and internal robustness in the plane in Sections [3] and |4] respectively. Subse- 
quently, we discuss internal robustness in the 3-dimensional space in Section [5] and 
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explore the robustness of some selected equilibrium classes in Section |6] Finally, in 
Section [7] we make additional remarks about smoothness and structural stability 
and use the previous results to formulate conjectures about the full (downward) 
robustness of equilibrium classes in 3 dimensions. 

2. Basic notations 

In this paper, we deal with convex bodies in the Euclidean spaces and M'', 
of dimension 2 or 3, respectively, where, by a convex body, we mean a compact 
convex set with nonempty interior. For a point p, we let \p\ be the Euclidean norm 
of p, and denote the Euclidean unit ball of the ambient space by B. 

We distinguish three subclasses of convex bodies: by Vn, On and /C„, we denote 
the families of n-dimensional convex polytopes, convex bodies with smooth {C°°- 
class) boundary, and convex bodies with piecewise smooth boundary, respectively. 

Let K E IC2 he a. convex body, and p £ intK. We say that q E bdiiT is 
an equilibrium point of K with respect to p, if the line passing through q and 
perpendicular to q — p, supports K. Clearly, if hdK is smooth at q, then this 
condition is equivalent to saying that g is a critical point of the Euclidean distance 
function z H- |z — p|, z £hdK. 

We call the equilibrium at q nondegenerate, if one of the following holds: 

• if bdif is smooth at g, then the second derivative of z t-^. |z — p|, z e bdif 
at q is not zero, 

• lihAK is not smooth at then both angles between p — q and one of the 
two one-sided tangent half lines of bd if at q are acute. 

Note that ii K € O2 ov K € 1^2, then this definition reduces to the usual concept 
of nondegeneracy in these classes. In the case of a smooth point, we call the 
nondegenerate equilibrium point q stable or unstable, if the second derivative at q 
is positive or negative, respectively. In the nonsmooth case, we call the equilibrium 
point unstable. The fact that the numbers of the stable and unstable equilibrium 
points of any K e /C2 are equal if K has only nondegenerate equilibrium points, 
follows from the Poincare-Hopf Theorem (cf. 2J). These two types of points form 
an alternating sequence in bd K. 

These definitions can be naturally adapted to convex bodies in O3 using the Eu- 
clidean distance function, (distinguishing three types of nondegenerate equilibrium 
points: unstable, saddle and stable points; depending on the number of negative 
eigenvalues of the Hessian) and also to convex polytopes in M.^ . We note that in the 
latter case, unstable, saddle and stable points are vertices, relative interior points 
of edges and of faces, respectively. For the definition of nondegeneracy in the piece- 
wise smooth case in higher dimensions, the reader is referred to :21j. If K has only 
nondegenerate equilibrium points, then the Poincare-Hopf Theorem yields that 

S-H + U = 2, 

where S', H and U denote the numbers of the stable, saddle and unstable points of 
the body. 

Throughout the paper, we deal with bodies that have only nondegenerate equi- 
librium points. 
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For simplicity, we denote the family of plane convex bodies K (z IC2, with S 
stable points with respect to their centers of gravity, by {S}. For convex bodies in 
O3 or in Vs, we may define the class {S,U} similarly, where S and U denote the 
numbers of the stable and the unstable points with respect to the center of gravity 
of the body. 

If K has S stable points with respect to some p G intK, let J^^{K,p) be the 
family of (convex) subsets of K such that any K' S {K, p) has strictly less than 
S stable points with respect to p. We obtain a similar notion, which we denote by 
J-<c{K), if the reference point is not fixed, but in the cases of both K and K' it is 
the center of gravity of the corresponding body. We define the (full) robustness of 
a planar convex body K S /C2 as follows. 

Definition 1. Let K G {5}. Then we define the downward robustness (or simply 
robustness) of K as the quantity 

min{area(ii' \ K') : K' G J'<{K)} 



p{K) 



area(-fi') 



In M'^, we may define p{K) in an analogous way. For brevity, we set ps = sup{p(_fi') : 
K e {S}} and ps,u = sup{p(if ) : K € {S, U}}. 

In Section [3] our aim is to investigate robustness with the reference point fixed, 
which we define below. Note that in this definition the reference point need not be 
the center of gravity of the body. 

Definition 2. Let K G K.2 and p G int K. Assume that K has S stable points with 
respect to p. We define the downward external robustness of K (or simply external 
robustness) with respect to p as the quantity 

, min{area(X \ K') : K' G T<{K,p)} 
P = ■ 

For simplicity, we set 

pf ^su^{p--{K,G):Ke{S}}, 

where G denotes the center of gravity of K. In M'^, we may define p'^^{K,p) and 
P's'u analogous way. 

In Sections |4] and [5j we examine internal robustness, which we define below. 

Definition 3. Let K G IC2 and p G int if. Assume that K has S stable points 
with respect to p. Let R{K, p) C denote the set of the points such that K has 
S stable points with respect to any point of R{K,p). The internal robustness of K 
with respect to p is 

tn(T^ ^ min{|g-p| : q j R{K,p)} 
perim K 

where perim if is the perimeter of K. 

Remark 1. By compactness, it is easy to see that if K has S stable points with 
respect to p, then p has a neighborhood U such that with respect to any z G U , K 
has S stable points with respect to z. Thus, if K has only nondegenerate equilibrium 
points with respect to p, then p^"-{K,p) > 0. 
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Remark 2. Definition [3] does not distinguish between upward and downward ro- 
bustness. If tlie latter are understood as tfie minimal distance necessary to move p 
to achieve increase/decrease of N then our definition refers to their minimum. 

Similarly like for full and external robustness, we set 

p's"^sup{r{K,G):Ke{S}}, 

wher G is the center of gravity of K. In R^, we define p™{K,p) and p^gjj similarly, 
by replacing perim K by the square root of the surface area surf K of the body. 

We note that it is easy to show that for each S and U the infimum of p{K,p), 
p'^^{K,p) and p^^{K,p) is zero in any class {S} and {S, U}. 



3. External robustness 



Our main result in this section is the following. 

Theorem 1. Let K (z IC2 contain the origin in its interior, and assume that K 
has S > 3 stable points with respect to o. Then 

with equality if, and only if K is a regular S-gon and o is its center. 

In light of the notations introduced in the previous section, we may reformulate 
this statement in the following way, for the special case that o is the center of 
gravity of K . 

Corollary 1. For any > 3, we have 



Ps 



ex „ tan I - g 



S'tanf ' 



and the planar convex bodies in {5*} with maximal external robustness are the 
regular iS-gons. 

First, we prove Lemma [T] 

Lemma 1. Let Q be a convex quadrangle with vertices a, b, c, d in counterclockwise 
order. Let the angles of Q at b and d be right angles, and let a G (0, tt) denote the 
angle at a. Set r = max{|fo — a|, jd — a|}, and for any plane convex body C ^ Q 
containing [a,b] U [a,d], let X{Q,C) = {z <E C : \z ~ a\ > r}. Then, the quantity 
p{Q,C) = ^ 

°arca(C'')'^'''' maximal over the convex quadrilaterals Q, with given 
angles, and over plane convex bodies C Q Q containing [a, b] U [a, d] if, and only if, 
C = Q and Q is symmetric about [a, c] . 

Proof. It is easy to see that (if we permit Q to be a degenerate quadrangle) the 
maximum of p{Q, C) is attained for some Q and C, hence, we show that if |6 — a| 7^ 
Id — a| or C 7^ Q, then p{Q, C) is not maximal. 

Without loss of generality, let r = |6 — a| > |c?— a|, and assume that Xc 7^ 0. 
Let z be the point of C with \z ~ a\ — r and farthest from 6, and let H denote the 
closed half plane containing [d, z] on its boundary, and a in its interior. 
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If [d, z] ^ bdC, then we may remove a part of C \ Xc, and thus increase /i. 
If Xc 7^ (-'^Q n H) then we may increase a,Tea{X {Q , C)) , and thus, increase fj,. 
Consequently, it suffices to examine the case that [d, z] C bdC and X{Q,C) = 
XiQ,Q)nH. 

Consider the case that C 7^ Q, which yields that z ^ [c,d]. Let L be the line 
passing through z and parallel to [c, d] , and H' be the closed half plane bounded 
by L and containing b. Let c' and d' denote the intersection point of L with [b, c] 
and [d,a] U [a, 6], respectively. Set Q' = conv{a, &, c', d'} if d' g [d,a], and Q' = 
conv{d',&, c'} if d' e [a,b] (cf. Figure[l]). Furthermore, set C" = Q' n if'. Observe 
that area(X (Q',C")) > area(X(Q,C)) and area(C" \ X(Q',C")) < area(C\ Xc). 
Thus, ^J.{Q, C) < n{Q', C), which yields that if /^(Q, C) is maximal, then C — Q. 




Figure 1 . An illustration for the proof of Lemma [T] 

Finally, we examine the case that C = Q and |6 — a| > |d — a|. Let [a',d'] be 
defined by the properties 

• a' e [a, b] and d' G [c, d] , 

• [a',d'] and [a,d] are parallel, and 

• \d' -a'\ = \b-a'\. 

Let Q' ^ C = conv{a\b,c,d'}. Observe that X{Q,C) C X{Q',C') and Q' C Q. 
Hence, C) is not maximal, and the assertion readily follows. □ 

Proof of Theorem^ Let the stable points of if be si, S2, . . . , ss in counterclockwise 
order in bd K. Let Ki be the closed set of points z oi K such that the position 
vectors Si, z and Si+i are in counterclockwise order around o. For any i — 1, 2, . . . , n, 
let r, = max{|si|, |si+i|} and set Xi = {z G Ki : \z\ > 7-^}. 

We observe that if if' C if has less than S stable points with respect to o, then 
for some value of i, K' O Xi = %. Indeed, any such K' has at most S — 1 stable 
points. On the other hand, if there is some point Zi € K' C\ X^ for every value of i, 
then the sequence \zi \ — |si|, |s2| — jzij, |z2| ~ IS2I, • ■ • , ~ \zs\ alternates, which 
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yields that the Euchdean distance function, defined on the points of bd C, has at 
least S local minima; a contradiction. Thus, noting that o ^ Xi for any value of i, 
we obtain that 

ej,, , min{area(Xj) : ?; = 1, 2, . . . , S"} 

P (^'0) = — . 

area(A) 

We intend to maximize this quantity over K G /C2- 

First, observe that if area(Xi) is not minimal for some value of z, then we may 
truncate Xi and thus decrease area(ii'). Hence, we may assume that aiea{Xi) is 
the same quantity for every value of i. Now, let = Z(si, o, Si+i). 

For each i with < tt, let K'^ — conv{xi, o, y^, Wi} be the right kite such that 

• Xi and Hi are on the half lines, starting at o, that contain Si and Si-)_i, 
respectively, 

• area(X-) = area{Xi), where X- — {z £ K[ : \z\ > \si\}. 

Note that by Lemma [l] we have aiea,{Ki) > area{K'^) for every z, with equality if, 
and only if Ki = K[. 

Case 1, if ai < it for every value of i. We let K' = K'^, x = area(X-), and 
note that K' is not necessarily convex, and that area{X^) is independent of i. 

Observe that p'^^{K,o) < ^^j^T^p^Ty, with equality if, and only if K = K'. Thus, 
it suffices to show that area(j<'') maximal if, and only if K' is a regular S'-gon, 
or in other words, if = ^ for every i. Since the maximum of ^j-ca^g') clearly 
attained, we show only that if ai ^ ai+i for some i, then this quantity is not 
maximal. 

Assume that oii ^ a^+i for some i. Set = \xi\ and r^+i ~ \xi^\\. Then 
2/ ai ai\ 2 (, "i+i 

*^"y-yJ^"'+i ^J' 



X — r. 

and 



area(ii'i) + area(i4ri+i) ~ tan + r^_^_l tan 

Using the notation a — ^ and /3 ~ the first condition can be transformed 

into the form 

tan a — a 



tan/3-^' 
which yields that 

area(is:,) + area(i^i+i) _ lZ"pZ$ tan (3 + tana ^ tan a tan /3 



A: 



tan a — a tan a — a tan /3 — /3 



We need only show that under the constraint that a +(3 is constant, A is minimal 
if and only ii a ~ B. But this indeed holds, since for the function f(a) — , 

f ' J \ / tan a — Q ' 

we have 

,, 2tana(Q;^ — sin^ a) 

/ (a) = 5—7; ^ > 

cos^ a(tan a — o;)'^ 

for every a e (O, f ), which yields that / is strictly convex. 

Case 2, if ai > tt for some value of i. Observe that in this case > tt for exactly 
one value of i. Let this value be S, and set J2i=i cti — < tt. Let K' — Ks U 
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(Uf=/^l)> and X = area(XO for some i ^ S. Then p^^{K,o) < ^j^^^. Using 
the argument of Case 1, we have that if a^j-ea^g') maximal, then K[, K2, . . . , K'g_^ 
are congruent right kites, with their angles at o equal to i^^z^i- According to our 
consideration, we have area(X5) — area(Xj') and area(/'ir5) > ai-ea{K'^) for any 
i S. Thus, in this case we have 

where P is a regular S'-gon, with the origin as its center. □ 

From the proof of Theorem [T] one can easily deduce Corollary [2j 

Corollary 2. Assume that K G IC2 has S stable points with respect to p € intii". 
Then there exists e — s{K) such that ii K' C K has less than S stable points with 
respect to p G intif', and area(iir\ K') — p'^^{K,p) area(i^) < e, then p''^{K',p) > 



4. Internal robustness in the plane 
Our main theorem is the following: 

Theorem 2. For any K G /C2 o,nd p € intii', if K has S > 3 stable points with 
respect to p, then p'^'^{K,p) < 1^ , with equality if, and only if, K is a regular S-gon, 
and p is its center. 

Similarly like in Section [sj we may reformulate Theorem [2] in terms of p™. 
Corollary 3. For any > 3, we have 

in _ J_ 

and the plane convex bodies K £ {S} with maximal internal robustness with respect 
to their centers of gravity are the regular S'-gons. 

In the proof we use Definition |4] and Lemmas [2] and |3] 

Definition 4. Let be a family of closed segments in M^, and ^ C a set. If 
for every [a, b] £ F there is a closed infinite strip, containing A and bounded by a 
pair of parallel lines La and Lb such that a € La and b £ Lb, then we say that A 
admits a strip cover by the elements of T . 

Lemma 2. Let P C 6e an S-gon of unit perimeter. If B = q + pB admits 
a strip cover by the sides of P, then p < ^ with equality if, and only if, P is a 
regular S-gon, and q is its center. 

Lemma 3. Let K £ IC2 have S > S stable points with respect to p £ intiC. Let 
P denote the convex hull of the unstable points of K . Then B = p -\- p^"{K,p)'B 
admits a strip cover by the sides of P. 



First, we prove Theorem [2] and then the two lemmas. 
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Proof of Theorem^ Let K G IC2 be of unit perimeter and let p G intif. Assume 
that K has S > 3 stable points with respect to p. Let P denote the convex hull 
of the unstable points of K. Since P <Z K, perimP < perimif = 1, with equality 
if, and only if, P = K . Thus, applying Lemmas [2] and [3j we immediately have 
p^"-{K,p) < On the other hand, if p"^{K,p) = then perimP = 1, and by 
Lemma [2] P is a regular S'-gon, and p is its center. From this, it readily follows 
that P = i^T is a regular S-gon, and p is its center. □ 

Proof of Lemma^ Observe that, for any side [a, b] of P, the width of any infinite 
strip, bounded by the lines La and L}, with a € La and b G P^, is at most |6 — a|, and 
here we have equality if, and only if La and Pf, are perpendicular to [a, b]. Thus, 2p 
is not greater than the length of a shortest side of P. Since perimP = 1, it readily 
implies the inequality p < jg. 

Assume that P — Then P is an equilateral S-gon, and for any side [a, 6] of P, 
the strip conv(La, Pf,), containing B, is perpendicular to [a, 6] and is circumscribed 
about B. Thus, the center q of P is on the bisector of the side. Since this holds for 
the bisector of each side of P, it follows that q is the center of the circle circumscribed 
about P. Note that if each vertex of an equilateral polygon lies on the same circle, 
then the polygon is regular, which immediately implies the assertion. □ 

Proof of Lemma^ Let K G IC2 be of unit diameter with S > 2 stable points with 
respect to the origin 0. Recall that the caustic of a smooth curve (also called evolute, 
cf. [7]) is the locus of the centers of curvature of the curve. It is well-known (cf. 
[19j). that a plane convex body with smooth boundary has degenerate equilibria 
only with respect to a point of its caustic, and that otherwise the number 2S of 
the equilibria of K changes (and in this case it changes by two) if, and only if, 
the reference point transversally crosses the caustic (in 3-dimensional space one of 
the two caustics). We use this observation in the following, more general form: 
p*" (K, 0) is the largest number p > such that 

(1) pintB contains no center of curvature at any smooth point of hdK, 

(2) pintB contains no point of any one-sided inner normal half line at any 
nonsmooth point of hdK. 

We parametrize bd if as i 1-^ where t G [0, 1]. For i = 1, 2, . . . , S, we denote 
the stable and unstable points of K hy Ui ^ Tit'i) a-nd Si — r{ti), labelling them in 
such a way that < ti < t[ < . . . < tn < t[^ < 1. Let P = convjiti, -U2, ■ • ■ , "s}- 

Clearly, for every i, 

(1) the tangent lines of hdK at Ui = r(ti) and at Si = r^ti) are perpendicular 
to the position vectors Ui and Si, respectively if Ui is a smooth point of 
hdK (note that Si is a smooth point by definition), and 

(2) the angle between Ui and any of the two one-sided tangent lines is acute, if 
Ui is not a smooth point. 

Thus, the angles of the triangle conv{o, at Ui and u^+i are acute, from 

which it follows that each side of P contains a stable point with respect to o. 

Observe also that as hdK is piecewise smooth and K has only nondegenerate 
equilibria, {r{t),r{t)) < at any smooth point with t G {ti,t[) (that is, \r{t)\ 
strictly increases), and it is positive at any smooth point with t G (t^,ti+i) (|r(t)| 
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decreases). This observation holds also for nonsmooth points, if we replace r{t) by 
any of the two one-sided derivatives of r{t). 

Let ai = /l{ui, o, Si), Pi = Z(si,o,Ui+i) and U = — Ui\. We distinguish 

three cases. 

Case 1, ai < ^ and /3i < ^. 

For simplicity, we imagine the segment [ui, Wi+i] as 'horizontal', and assume that 
o is 'below' this segment (cf. Figure [2]). Let L and L' denote the two lines, parallel 
to [o. Si], that pass through the points Ui and u^+i, respectively. Let B and B' 
be the closed disks, with o as their centers, that touch L and L', respectively. We 
show that both disks contain a center of curvature, or intersect a one-sided inner 
normal half line at a nonsmooth point of the arc ^([ti, ii+i])- This clearly implies 
that p*"(_fi', o)B is contained in a strip between two parallel lines passing through 
Ui and Wi+i, respectively. 




Figure 2. The region in Case 1 of the proof of Lemma |3] 

We show only that B contains a point of the caustic or a one-sided inner normal 
half line, as for B' we may repeat the same argument. 

First, note that if Ui is not a smooth point of hdK, then the right-sided inner 
normal half line at Ui intersects B. Thus, we may assume that Ui is a smooth point 
of bdi^:. 

Let w denote the tangent point of B on L. Observe that since < ^, u,; is 
upwards from w. Let D denote the region, containing r^, that is bounded by the 
union of the following arcs: 

• the closed half line in L, emanating from Ui and containing w; 

• the points r{t) with t e [^4,^^]; 

• the segment [o, Si] ; 
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• the closed half ray in the line containing [o, Ui] , starting at o and not con- 
taining Ui- 

This region is shown as a dotted domain in Figure [2] 

Let i?t denote the inner normal half line of — r{[ti, t'^\) at r(t), if it exists. Since 
{r{t),r{t)) < for every t e {ti,t'^), we have that in this interval Rt n [o,Si] = 0. 
Thus, Rt C D for every t € [ti, t ■]. Thus, the assertion readily follows if gi contains 
a nonsmooth point. 

Now, consider the case that gi is a smooth curve, and let and denote the 
centers of curvature at Ui and Si, respectively. Observe that qi e relint[wi,o] and 
o G relint [si, fi]. Furthermore, we have qi S B, or S B, or that g^ and 
are not in the same connected component of D \ B. Since in the first two cases 
the assertion readily follows, we may assume that qi and are not in the same 
connected component. 

If the curvature of bd C is not zero at any point of gi , then, as in this case the 
caustic of gi is a continuous curve in D that connects qi and r^, B contains a center 
of curvature of gi. Consider the case that the curvature of gi is zero at some point. 
Let t g [tijt'i] be the smallest value such that the curvature of gi is zero at r{t). 
Then, the caustic of r{[ti, i)) is a continuous curve in D that connects Ui to a point 
in the unbounded component of D \ B, and thus, contains a point in B. This 
finishes the proof in Case 1. 

Case 2, ai + Pi < tt, but one of the two angles is at least |. Let, say, ai > ^ (cf. 
Figure [3]), which readily implies that /3i < ^. Let d denote the distance of o and the 
line through u^+i and parallel to [o, Si]. This distance is the length of the segment 
[ui+i, y], where y is the orthogonal projection of u^+i on the line containing [o, Si]. 
Let X be the intersection point of and [o, s^]. 




Figure 3. An illustration for Case 2 of the proof of Lemma [3] 

Observe that by the argument in Case 1, we have that d > p^"^{K,o). On the 
other hand, p"^{K,o) < \ui\. Indeed, if Ui is a smooth point of hdK, then |ui|B 
contains the center of curvature at Ui, and if Ui is not a smooth point, then it 
intersects the right-hand side inner normal half line at Ui . 

Consider the case that — a;| < a;| (cf. Figure[3|. Then, since (i^, o) < 
d, choosing L and L' parallel to [o,Si] and satisfying Ui G L, Ui+i € L' , we have 
that p^'^iK, o)B C conv(L U L'). 
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Now we examine the case that \ui^i — x\ > ju^ — Let L and L' be the hnes 
perpendicular to Uj such that Ui G L and Ui+i € L' . We show that p^"{K,o)'B C 
conv(L U L'). Let z be the intersection point of [ui, u^+i] with the hne containing 
o and parahel to L. Observe that since > f, we have z G Thus, 
— z| > I Mi — z I, which yields that |ui|B C conv(iUi'). Since \ui\ > p'^^{K,o), 
it readily implies that p^"'{K,o) C conv(L U L'). 

Case 3, ai + f3i > tt. From this condition, (using the terminology of Case 1) it 
immediately follows that o is 'not below' the segment [ui, Ui+i], and thus, o ^ int P 
(cf. Figure |4]). Observe that the condition of Case 3 may hold for at most one value 
of i. According to the previous cases, we have that for every j ^ i, p™{K,o)'B is 
contained in a strip bounded by two parallel lines passing through Uj and Uj+i, 
respectively. We may even observe that p is contained in the halves of these strips, 
bounded by the corresponding sides of P, that overlap P. 




Figure 4. An illustration for Case 3 of the proof of Lemma [3] 

We show that there is such a strip for [ui,Ui4.i] as well. Indeed, let Hi_i denote 
the half of the infinite strip, belonging to [ui-i^Ui], that is bounded by [ui-i,Ui\ 
and contains o. Define Hi^i similarly for [m^+i, Mi-|_2]. Let Li denote the ray in 
bd-ff; starting at Ui, and L^+i denote the ray in bd-ffi+i starting at u^+i. 

If Ui ^ -ffj+i, then the strip, bounded by L^+i and its translate starting at Ui, 
contains -ffi+i, and thus also p"'^ {K , o)'B^ . If Ui+i ^ Hi, we may apply a similar 
argument. Thus, we may assume that [ui,Ui+i] C HidHi^i (cf. Figure |4]). Let cjji 
and (pi+i denote the angle between Li and [ui, u^+i], and the angle between L^+i and 
[ui, Ui+i], respectively. Note that since u^+i £ Hi, we have (j>i < Z(ui_i, Ui+i, Ui) < 
Z(iti+2j Mi) < ^. We may obtain similarly that 4'i+i < ^. Hence, the strip 
bounded by the lines perpendicular to [ui,Ui+i] and passing through its endpoints 
contains Hi n Hi+i, and thus, also p™(i4r, o)B^. □ 

5. Internal robustness for 3-dimensional convex bodies 

In this section we partly generalize the results of Section [4] for convex polyhedra. 
Our main result is as follows: 
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Theorem 3. Let P he a regular polyhedron with S faces, U vertices and H = 
S + U — 2 edges, and let a be the center of P. Let P' be a convex polyhedron with 
S faces, U vertices and H edges, each containing an equilibrium point with respect 
to some q G int P' . Then 

p'-{P',q)<nP,o), 
with equality if, and only if P' is a similar copy of P , with q as its center. 

Remark 3. Clearly, we do not need to require formally that the numbers of faces, 
edges and vertices, stable, saddle and unstable points of P' are equal to the corre- 
sponding quantities of P. We leave it to the reader to show that if, say, P' has S 
stable points, U unstable points and H edges, then it has S faces, U vertices and 
H saddle points. 

In the proof we use the following theorem of Dowker (cf. jl5j). 

Theorem 4 (Dowker). Let C be a unit circle, and G be an n-gon circumscribed 
about (equivalently, containing) C. Then: 

• area(G) is minimal if, and only if, G is a regular n-gon circumscribed about 
C, 

• denoting by a„ the area of a regular n-gon circumscribed about C , the se- 
quence {a„} is strictly convex; namely, for any n~2 > k > 0, On-k+O'n+k > 
2a„. 

Proof of Theorem^ For simplicity, assume that surf(P') — surf(P) = 1, and let 
R denote the radius of a circle inscribed in a face of P. Furthermore, for any face 
F oi P' , let rp denote the inradius of F, and let T denote the family of the faces 
of P'. Clearly, p'=^(P', q) < r ^ mm{rF : F E T}. 

By Theorem |4j we have that r < R, with equality if, and only if, all the faces 
of P contain a stable point, and the faces of P are congruent regular polygons. 
We note that since replacing an n-gon by a regular n-gon, and replacing a regular 
(n — fc)-gon and a regular (n -I- fc)-gon by two regular n-gons does not change the 
total number of edges, we have that, in the case that r = R, the faces of P' are 
congruent to the faces of P as well. 

On the other hand, p"\P' , q) — r implies also that all the lines, orthogonal to a 
face and passing through its incenter, meet at q. Furthermore, if the faces F and 
F' are joined by an edge, then, clearly, the distances of q from the centers of F 
and F' are equal. Thus, q is at the same distance from any face, which means that 
the face angles between any two faces joined by an edge are equal. From this, we 
obtain that the vertex figures of P are congruent, which yields that P' is a regular 
polyhedron. Hence, the assertion readily follows. □ 

The following example shows that Theorem |3] is false without the condition that 
the numbers of the edges of P and P' are equal. 

Example. Let P be a regular tetrahedron of unit surface area with center o. 
Truncate P near a vertex, in such a way that does not change the numbers of the 
three types of equilibria of P, and the truncated part does not intersect the incircle 
of any face of P, and denote the truncated polyhedron by P'. Then P' has the 
same numbers of stable, saddle and unstable points with respect to any point of 
int(p™(P,o)B), but surf(P') < surf(P) = 1. Thus, /"(P',o) > p™(P,o). 
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6. Full robustness, and its variants, in 3 dimensions 

In this section, first, we show that the full robustness of equilibrium classes 
{S, U} with S,U < 3 is maximal, and then we apply our method for some other 
types of robustness. 

Theorem 5. We have pi2 = P21 = P22 = 1- 

We start by introducing a definition and proving some lemmas. 

Definition 5. For a convex body iv' in K^, a hounding box of -fC is a brick circum- 
scribed about K. 

Lemma 4. If for some hounding box of a convex body K cM.^ with edge lengths 

a < b < c 

• we have 66 < c, then K has at least two unstable points with respect to its 
center of gravity, and if 

• we have 3a < b, then K has at least two stable points. 

Proof. Let B be a bounding box of K with edge lengths a<b <c. 

First, we show that the distance of the center of gravity of K from any face of 
K is at least one quarter of the width of B in that direction. Consider two parallel 
faces of -B, say those at the distance c from each other. For simplicity, we assume 
that the center of gravity of K is the origin o, and these two faces are in the planes 

z = Zi > 0, and z = —Zi < for some 21, 2:2 G K with zi + Z2 = c. 

Let Kq be the intersection of K with the {x, y)-plane and let pi be a point of 
K in the plane z = zi. Let C be the infinite cone with p as its apex and Kq as 

its base, and let K' be the intersection of C with the infinite strip bounded by the 
planes z = Zi and z = —zi- Then, by the definition of the center of gravity of K, 
we have 



= / z dx dy dz> z dx dy dz 

J(x,y,z)eK J (x,y,z)eK' 



3 

Zi - -c. 



From this Z2 > f (and similarly zi > |) readily follows. 

Now, assume that for the edge lengths of B we have 6a < 66 < c. Observe that 
for any q e relbdii'o, we have \q\ < | \/a^ + 6^ < -^c. On the other hand, for any 

point p G K in z = Zi OT in z = —Z2, we have > |. Since | > ^c, it means that 
the Euclidean distance function has at least two local maxima, one with a positive 
and one with a negative 2;-coordinate. 

If 3a < 6, then we may apply a similar argument. □ 

Lemma 5. There is some /x > such that for any A e M, there are convex bodies 
Ki2{X) € {1,2}, /C2i(A) G {2,1} and /C22 (A) € {2,2}, and their bounding boxes 
Bij{X) with edge lengths aij{X) < bij{X) < Cjj(A), where ij € {12,21,22}, such that 

' '^<- (,i2(A) vol(Bi2(A)) ^ ^' 
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Class {1,2} 



Class { 2,1: 




2r 


d 


2r 









Figure 5. Examples of convex bodies in classes {1, 2} and {2, 1}. 



Proof. Figure [5] shows two convex bodies in {1, 2} and {2, 1} with bounding boxes of 
edge lengths a — b = 2r and c = d + ir and a = 2r and c> b > d, respectively. The 
body in {1, 2} is a suitably truncated cylinder. We construct the body in {2, 1} by 
tilting two, originally coincident circular discs by a small angle. The section of the 
two tilted discs with the plane containing their centers are two straight line segments 
and we connect those by tangentially attached circular arcs. We remark that the 
convex hull of a Dupin 'needle' Cyclide appears to be in the same equilibrium class. 

Let -^"12 (A) and ^"21 (A) be these bodies, respectively, with some d > 2rA. Then 
the edge lengths of the given bounding boxes of the bodies satisfy the conditions 
of the lemma. To obtain the bodies Ki2{X) and K2i{X) with smooth boundaries, 
we may apply the standard smoothing algorithm, described, for instance, in [TU]. 
Then it is an elementary computation to show that the limits of vo1(J3i2(a)) 
voi(K2i(A)) strictly greater than zero. 



and 



vol(B2i(A)) 

As K2i{X), we may simply select an ellipsoid with semi-axes a 
-4A2. 



1, 6 = 2A, 
□ 



Now we return to proving Theorem [5] 

Proof. We start with the obvious remark that if any equilibrium class {S,U} con- 
tains, for every e > 0, a convex body with downward robustness at least 1 — e, then 
ps^u — 1- Hence, using the notation of Lemnia[5j we need only show that 

lim p(ifi2(A)) - lim p(/^2i(A)) = lim p(X22(A)) = 1. 

A— foo A— >-oo A— >-oo 



We start with the class {1,2}. The only way to reduce the number of the equi- 
librium points of _ftri2(A) is to truncate the body in such a way that the remaining 
object is in class {1, 1}. For brevity, we let 7^ii(A) denote a convex body obtained 
in this way by a truncation of "almost" minimal relative volume (note that a trun- 
cation of minimal relative volume may yield a body with degenerate equilibrium 
points). Observe that at least one bounding box of Kii(X) fits inside the bounding 
box of Ki2{X) described in Lemmajij Thus, by Lemmajijand as vo!(ai2(A)) — ^ ^'^^ 
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every value of A, we have 

lim p(ifi,(A)) = l- lim "°!i'!"!'ii =l. 

A^oo A^oo VOl(iiri2(A)) 

The argument for the class {2, 1} runs in an analogous manner. 

In case of the class {2,2}, there are three alternative ways to produce a convex 
body with less than six equilibrium points: we may obtain a truncated body in one 
of the classes {1,2}, {2,1} or {2,2}. By Lemma [3} for any bounding box of the 
truncation with edge lengths a < 6 < c, we have 66 > c or 3a > b. Thus, applying 
an argument similar to the one in the previous cases, we have limA->.oo p(-^22(A)) = 
1. □ 



Until this point we treated stable and unstable points in a similar manner, ro- 
bustness corresponded to truncation resulting in the reduction of the number N 
of equilibria, regardless of their type. It is natural, and, as we will see, also useful 
to ask for partial robustness, i.e. the (relative) volume of a truncation necessary 
to reduce either S or U, (the numbers of stable and unstable points, respectively,) 
resulting in the notion of S"- robustness and L/- robustness, denoted by p^, p", respec- 
tively. Naturally, we have 

p = min{p^p"}. 

We can define the S- and robustness of equilibrium classes, denoted by pl^^p^j, 
respectively, in a natural way, and we can immediately see that 

because it is not possible to further reduce the number of stable and unstable 
equilibria in these classes. Beyond this trivial comment, with ideas very similar 
to the ones in the proof of Theorem [5j it is easy to obtain one additional result 
about the partial robustness of two infinite families of equilibrium classes which we 
formulate in the following theorem. 

Theorem 6. If n > 2 then n = 2 — 1- 

Proof. Let and if„2 be a regular, planar n-gon, and a straight, infinite prism 
with a regular n-gon as its base. We leave it to the reader to show for every A > 
the existence of the convex bodies K2n{^) & {2, n} and Kn2{X) € {n, 2} with some 
bounding boxes satisfying the conditions in the first two parts of Lemma [sj These 
convex bodies approach K2n and Kn2, respectively, as A — >■ oo. Then, to prove the 
assertion, we may apply Lemma [4] and follow the idea of Theorem [5] □ 

Theorem |6] leads to the following Corollary. 
Corollary 4. If ti > 2 then p2,« = P2,n Pn,2 = pfi,2- 



7. Remarks and open problems 



In this paper we investigated the robustness of a convex solid with TV equilibrium 
points. Beyond the original definition of 'full' robustness we also defined internal 
and external, in the latter case also upward and downward robustness. 
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7.1. Smoothness and structural stability. Before formulating conjectures and 
raising questions we make some remarks about smoothness and structural stability. 

Remark 4. We have examined convex bodies with piecewise C°°-class bound- 
aries. Nevertheless, all our results (and proofs) can be applied also for bodies with 
piecewise C^-class boundaries. 

Remark 5. Let > 3. For every e > 0, there is a (5 = S{e, S) > 0, with 
lime_>.o+o (5 = 0, such that if i^T G /C2 has n stable points with respect to p G int K, 
and p'^"'{K,p) > 5^ — e, then the Hausdorff distance of K and a regular S'-gon, with 
p as its centre, is less than S. 

Remark 6. Remark [s] holds also if we replace p^"{K,p) with p'^^{K,p). 

Remark 7. For every S > 3, the maximum of p"^{K,p) over IC2 can be approached 
by regions with smooth boundaries as well; or in other words, = sup{/9™(ii', p) : 
K £ 02,p E inti^}. To show this, it suffices to replace the boundary of a regular 
S'-gon near the vertices by suitable elliptic arcs. This provides a C^-class curve as 
the boundary of a convex region, which, after applying a smoothing algorithm like 
in [1^, yields a convex region K £ O2, with internal robustness "almost" equal to 

7.2. Conjectures and open questions. Since we know [TT] that the first nonempty 
class in the plane is {2}, trivially we have 

P2 = 1- 

Based on Theorems [T] and [2j we can formulate the next conjecture. 
Conjecture 1. p„ — ^(regular n-gon) if n > 2. 

In 3 dimensions the situation is less transparent, however, there are many inter- 
esting questions. As a modest generalization of Theorem[3]we propose the following. 

Conjecture 2. Theorem[3]is also valid for downward external robustness, and also 
for downward full robustness. 

A bolder generalization of Theorem [3] refers to equilibrium classes. We call the 
equilibrium classes containing platonic solids (classes {4, 4}, {6, 8}, {8, 6}, {20, 12} 
and {12, 20}) briefly platonic classes. An affirmative answer for Conjecture [2] would 
suggest that in the platonic classes platonic solids have maximal downward full 
robustness. This we pose as an open question: 

Problem 1. Prove or disprove that in the platonic classes platonic solids have 
maximal downward full robustness. 

So far we addressed the robustness of nine classes in M'^: in Section [g] we dis- 
cussed the four with indices less than 3 (Theorem [s]) , and here we formulated a 
conjecture related to the five platonic ones. In addition. Theorem [6] and Corol- 
lary [4] we addressed the S- and J7-robustness of two infinite class families; here we 
complement the latter results with 

Conjecture 3. P2,n = Pn,2 = Pn- 
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The underlying geometric idea is that in classes {2, n} and {n, 2} the very same 
solids (described in the proof of Theorem [6]) which have 'maximal' S- and U- 
robustness, also have 'maximal' overall robustness. In general, one would expect 
that a 'uniform distribution' of equilibrium points is associated with maximal ro- 
bustness. If all equilibrium points were of the same type, this would resemble the 
Tammes Problem (23j . however, since we have different types of equilibria, this 
analogy is incomplete. While our understanding of the robustness of general equi- 
librium classes is rather limited, the following conjecture is strongly suggested by 
Corollary [2] and is consistent with all our previous findings and conjectures: 

Conjecture 4. If i > fc and j > I then pij < pk.i- 

In general, it is an interesting question whether pij = pj^i. From the point of 
view of applications and numerical experiments it might be of interest to study 
partial robustness, measuring the difficulty of reducing separately either S or U. 
For example, based on the above arguments and conjectures, we expect in classes 
{2,n} the partial S'-robustness to be much smaller then the partial J7-robustness 
and the inverse statement applies for classes {n, 2}. 

Another 'reduced', however potentially important version of robustness can be 
defined by admitting only truncations by planes (or straight lines in 2 dimensions). 
This case is not only interesting because numerical experiments appear feasible 
but also this is the only truncation of a convex body where both resulting objects 
remain (weakly) convex. This naturally leads to the statistical study of the evolu- 
tion of numbers of equilibria in a population of convex solids which are generated 
by subsequent planar truncations and where both convex solids resulting from a 
trucation are considered. 

This idea also suggests an alternative definition of robustness. So far we consid- 
ered minimal truncations changing the number of equilibria, nevertheless it is also 
possible to regard the average magnitude of such truncations. While it is not clear 
how to define a measure directly on the space of general convex truncations, we 
can obtain a possible approximation by admitting successive truncations by (hy- 
per)planes. The space of the latter is a special type affine Grassmannian manifold 
with a natural measure fi, invariant under the motions of the embedding Euclidean 
space (cf. |20)). and can be extended in a natural way to a measure /i„ on the space 
of n successive planar truncations. We call a truncation nontrivial if it intersects 
K and we call it neutral, if it leaves the equilibrium class invariant. Then we can 
define the n-th order average robustness pn{K) associated with a convex solid K as 
the measure of the set of n, subsequent neutral truncations divided by the measure 
of the set of n, subsequent nontrivial truncations. Such a concept may be useful in 
computations and practical applications, also, it would be of interest to see whether 
solids with maximal average robustness also have maximal (downward) robustness 
in the original sense. 

One possible, practical choice for visualizing average robustness is to assume 
the above-mentioned natural measure on the space of truncations and then plot 
the percentage of truncations resulting in a given change AiV of the number of 
equilibrium points versus the relative volume corresponding to the truncation. This 
visualization admits the representation of an arbitrary convex solid of arbitrary 
dimension on the unit square. The area on the unit square corresponding to neutral 
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truncations (AiV = 0) could also serve as a definition of average robustness. In 
Figure |6] we show this plot for the square: the lower, middle and upper regions 
represent the truncations resulting in a polygon with S = 3, S ~ 4 and 5* = 5 
stable points, corresponding to AA^ = —1, 0, 1, respectively. 



Raliaof vwicalksns 
w'gh afnced nwnber 
of stable poii^s 




0.2 0.4 0.6 O.S 

Relative area of the trtmcatton 



Figure 6. Truncations of a unit square with one line. 

As pointed out in the introduction, our results may help to understand natural 
abrasion processes. Theorems [5] and [6] indicate that if we study the statistical 
outcome of a natural abrasion process, the fraction of solids in classes {1, n}, {n, 1} 
would be rather low as transition from neighbour classes {2, n}, {n, 2} is prohibited 
by the maximal S- and JJ-robustness of the latter. Indeed, field studies [T^] 
appear to confirm this: in coastal reagions the percentage of pebbles in classes 
{l,ri},{n, 1} was found to be below 0.1%. 
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